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Abstract 

We consider proton stability based on £"6 inspired extra U(l) model with S4 x Z2 flavor symmetry. 
In this model, a long life time of proton is realized by the flavor symmetry in several ways. We classify 
the suppression mechanisms of proton-decay and explain how the flavor symmetry works. There is an 
interesting solution, such as, in a special direction of vacuum expectation value (VEV), baryon number 
violating interactions are canceled. In the case that the suppression of proton decay is realized by the 
appropriate size of VEV, the allowed region of VEV exists when the exotic quark mass is O(TeV). From 
the constraint for the life time of exotic quark, the right handed neutrino mass should be in narrow range 
around lO^^GeV. 
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1 Introduction 



Supersymmetry is an elegant solution of hierarchy problem of Standard Model (SM) [1] and gives a new 
view point of generation structure of leptons and quarks. As a simple supersymmetric extension of SM 
suffers from non-conservation of baryon number, we must introduce R-parity symmetry in order to avoid too 
rapid proton decay. This means we can not construct any consistent superpotential based on only SM gauge 
symmetry. Even if we introduce R-parity, the superpotential of minimal supersymmetric standard model 
(MSSM) is not perfect one. The superpotential of MSSM suffers from //-problem such as we must tune 
the scale of //-parameter to be O(TeV), which is much smaller than Planck scale. Therefore the R-parity 
symmetry should be replaced by other symmetry. The information what symmetry should be introduced 
may be extracted from the structure of Yukawa interactions because these interactions are derived from 
superpotential. 

The appropriate start point is given by introducing an additional C/(l) gauge symmetry to forbid //-term 
[2]. In this frame work, several new superfields such as singlet S, exotic quarks G.C^, must be introduced 
to cancel gauge anomaly. Then the baryon number violating interactions in superpotential are replaced by 
single exotic quark interactions, which make it easy to suppress proton decay by the new symmetry. 

Considering the Yukawa interactions, we can guess about which symmetry we should introduce. Strangely, 
the mixing angle 6*23 of Maki-Nakagawa-Sakata (MNS) matrix is almost maximal. Many authors discussed 
non-Abelian discrete flavor symmetries to understand the structure of the MNS matrix. The flavor symmetry 
may be a good candidate for replacing R-parity symmetry. At previous work, we explained S'4 x Z2 flavor 
symmetry not only realizes maximal mixing angle ^23 but also suppresses proton decay based on SU{2i)c x 
SU{2)w X U{1)y X U{l)x X U{l)z gauge symmetry [3]. Therefore there is a deep connection between the 
proton stability and the existence of generation in supersymmetric model. 

In this paper we give more detailed estimation of proton life time and classify the mechanisms to suppress 
proton decay. At first we give definition of our model in section 2, and give classification of suppression 
mechanisms of proton decay in section 3. In section 4, we check that our superfield assignment realizes 
neutrino mass square differences and MNS matrix including new experimental value of ^13. We explain the 
origin of U{l)z breaking scale in section 5. Finally we give conclusion of our analysis in section 6. 

2 S4 X Z2 flavor symmetric extra U(l) model 

At first we explain the basic structure of our model. We extend the gauge symmetry to SU{3)c x SU{2)w x 
U{1)y^U{1)x^U{1)z which is the maximal subgroup of Eq. In order to cancel gauge anomaly, we must add 
new superfields, such as SM singlet S, exotic quark G, G"^ (hereafter we call them g-quark) and right handed 
neutrino (RHN) N". We can embed these superfields with MSSM superfields Q,U'', D", L, E", , into 
27 of Eq [4]. As the singet S develops VEV and breaks U{l)x gauge symmetry, O(TeV) scale /t-term is 
induced naturally. In order to break U{l)z and generate a large Majorana mass of RHN, we add SM singlet 
$, The gauge representations of superfields are given in Table 1 [3]. 
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Table 1: G2 assignment of superfields. Where x, y and z are charges of U{l)x, U {1)y and U{l)z, respectively. 
Y is hypercharge. After the gauge symmetry breaking of three U{l)s, R-parity symmetry R ~ exp[i7r(3a; — 
8y + 15z)/20] is unbroken. 

Under the gauge symmetry given in Tabic 1, the renormalizablc superpotential is given by 

W = Y'^H'^QU^ + Y°QD'^H^ + Y^H^LE" + Y^H^^LN^ + Y"^N''N^ + XSH^^H" + kSGG'' 
+ M'^'^" + yQQGQQ + Y^"G^UW^ + Y^'^GE'=U'= + Y^'^G''LQ + Y^°GN'D''. (1) 
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In this supcrpotcntial, unwanted terms are included in the second hue. The first term of the second hue is 
tire mass term of singlets $, 'I''^ which prevent singlets from developing VEVs. The other five terms of the 
second line are single g-quark interactions, which break baryon and lepton number and induce rapid proton 
decay. In the first line, we must take care of the fiavor changing neutral currents (FCNCs) induced by extra 
Higgs bosons [5]. Therefore the superpotential Eq.(l) is not consistent at the present stage. 

In order to stabilize proton, we introduce S4 x Z2 flavor symmetry. If we assign G, C to S4 triplet 
and quarks and leptons to doublet or singlet, the single g-quark interaction is forbidden. However, as the 
g-quark must never be stable from phenomenological reason, we assign to iS'4 triplet to break the flavor 
symmetry slightly. In order to realize the maximal mixing angle of 623 in the MNS matrix and suppress the 
Higgs-mediated FCNCs, we assign the superfields in our model as given in Table 2 [6]. 

In the non-renormalizable part of superpotential, the single g-quark interactions which contribute to the 
g-quark decay are given as follows 

yQQ yUD EU QL 

Mp Mp Mp Mp 

The detail of Wb depends on the Z2 charge assignment of pq and pg . 





Qi 


O2 


Q3 


U{ 








^2 


^3 


S4 


1 


1 


1 


1 


1 


1 


1 


1 


1 


Z2 


Pq 


V<i 


V<i 




Pq 


Pq 


Pq 


Pq 


Pq 




El 


El 


Et 


U 


L3 






Hi' 




S4 


1 


1 


1' 


2 


1 


2 


1 


2 


1 


Z2 


+ 




+ 






+ 






+ 










S3 


Ga 






$3 




S4 


2 


1 


2 


1 


3 


3 


2 


1 


3 


Z2 




+ 




+ 


Pa 


Pg 


+ 


+ 


+ 



Table 2: S4X Z2 assignment of superflelds (Where the index i of the S4 doublets runs i ~ 1,2, and the index 
a of the 5*4 triplets runs a = 1, 2, 3. The charges of quarks and g-quarks; Pq and Pg, take ±.) 



3 Classification of the suppression mechanism of proton decay 

Depending on the Pq and Pg, the flavor symmetry works in different ways to suppress proton decay. In this 
section we classify the suppression mechanism and estimate the allowed parameter range. 

(1) Leptoquark solution 

If we assign (pq,pg) = (+, — ), it results y*^*^ = y^^ = in Eq.(2). In this case we can assign lepton 
number (L) and baryon number {B) of G to {L,B) = (-1-1,1/3) and those of G'^ to {L,B) ~ (—1,-1/3), 
then the baryon number is conserved. Therefore proton becomes stable and VEVs of $, are not bounded 
from above. In this solution, G, G^ are well known as leptoquarks. Note that the VEVs of $, $^ are bounded 
from bellow, because the life time of g-quark must be shorter than 0.1 sec, otherwise the success of BBN is 
spoiled [7]. 

If we assign {pq.pg) = (—,—), then it results y^^ = y^^^ = y^^ = and g-quark becomes leptoquark 
which couples only to right handed charged leptons Ef,E^. Therefore the decay mode of our leptoquark 
depends on the flavor charge assignment. 

(2) Cancellation solution 

If we assign {pq,pg) = {+, +), it results y*^^ = 0. In this case, as we cannot deflne the lepton number and 
baryon number of g-quark, these numbers are not conserved. Here we investigate proton decay interactions 
which are induced by scalar g-quarks exchange. In the present flavor assignment, the superpotential which 
contributes to the proton decay is given by 

= -^C/,^i;3^[V3(G2$^-G3$^)$2-(G2*^ + G3$§-2Gi$^)$i] 

+ Y^QaQbM^G, + ^G2 + nG3) + ^KD^,^3mGl + $^G^ + <i>^G^) 
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+ ^QaQb[Vi{G2'^l - G3$§)$1 + (G2$^ + G3$§ - 2Gi$D*2] 
Mp 

+ ^C^:^fc[V3(G^$=2 - G^4>^)<I>i + (G^$^ + G§$§ - 2G^$5)$2]. 



(3) 



Note that the contribution from is omitted because it is shown to be r lepton in section 4 and does not 
contribute to proton decay. Integrating out scalar g-quarks, we get the effective four-Fermi interactions as 
foUows 
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where /i"^ = 63 and mean scalar g-quark mass Mq and dimcnsionless coefficients CgGtCgc a-re defined in 
appendix. Note that the masses of scalar g-quarks are degenerated due to the 6*4 symmetry. Interestingly, 



it results Xatc = A' 



= in the special VEV direction such as 



(5) 



In this case proton decay is forbidden. This means the contributions from three scalar g-quarks are canceled. 
(3) Suppression solution 

In the case that there is no remarkable cancellation, the size of ($) must be in appropriate region where 
the constraints for proton and g-quark life time are satisfied at the same time [8] . In order to suppress proton 
decay, the VEV of $ must not be too large. At first we estimate the upper bound of the VEV. As there 
are many unknown parameters in Wb and the VEV direction of $, is also unknown, we make several 
assumption for simplicity. At first, we assume there is no mixing between scalar g-quarks G and C and 
put Ggg = 1, Ggg<= = 0. Next, we change the assignment of Q, W^, D'^ and G, G'^ to {S4, Z2) = (1', +) and 
(54, Z2) = (3, -t-) respectively in Table 2 and replace the superpotial in Eq.(3) by 
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(6) 



where we assume ya = Zah = and the contribution from y^^ is omitted. Finally we tune the VEV direction 
as follows 
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Including the renormalization factor Aj^p, the effective four- Fermi interactions at IGeV is given by [9] 

^4 
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Arf = {Alp)siAI,p)siARF)L, (9) 
where we estimate the short distance part of A^p by the 1-loop renormalization group equations as follows 

, dlnZn 16 
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Here only QCD correction is accounted. This approximation is not bad because the beta function of the 
couphng constant of strong interaction vanishes at 1-loop level in our model, which makes the contribution 
of as dominant in the RGEs of Za and yah- Solving Eq.(lO) and Eq.(ll), if we put as{Mz) = 0.118, we get 



The long distance part is given by [10 

{Arf)l 
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(12) 
(13) 

(14) 
(15) 



For the case that the final state includes , the strongest experimental bound for the partial decay width 
of proton is given by p — > tt*^ + /i+ . For simplicity, we assume zi = Z2 = 23 = 1 for the mass eigenstates 
and tune ybc to normalize four-Fermi interaction as follows 



Ml 



From this Lagrangian, the proton decay width is given by [11] 
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where F and D are chiral Lagrangian parameters, is hadronic matrix element, is pion decay constant 
and and nip are masses of pion and proton. If we put 

F = 0.47, L» = 0.80, ap = -0.012 GeV^ = 130 MeV, m^o = 135 MeV, mp = 940 MeV[12], (18) 
then we get 



r(p^7r" + /i+) = (5.01 X lO-^^GeF) 
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From the experimental bound t{p 7r° + /i+) > 473 x 10^°[years] [13], the upper bound for VEV is estimated 
as follows 
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Next, we estimate the life time of g-quark under the assumption that g-quark is lighter than scalar g-quark. 
For simplicity, we assume the g-quark can decay only into smuon but not into stau or squarks. With this 
assumption, the g-quarks decay through the following interaction 



{A^rf)sV' 
Mil 



(u= + c^ + <^)/i=(5i+.g2 + .93), 



(21) 



from which one can see that g-quarks have the same life time. Requiring the life time of Qa is shorter than 
0.1 sec as follows 
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we get 

where Mg is g-quark mass. Hereafter we assume the approximation Mg = Mq is held for simpHcity. From 
Eq.(20) and Eq.(23), the allowed region for V is given by (see Fig.l) 

nOOOGcV\ ( ,,_27 ( Mg 



2.25 X 10"^" < < 2.96 x 10"^' . (24) 

V Mg ) \Mp) VlOOOGeVy ^ ' 

This inequality holds when the mass bound, 

Mg > 1.96 TcV, (25) 
is satisfied. For example, if we put Mg = 10 TcV, allowed region for V is given by 

Note that the factor of this constraint should not be taken seriously, because there is large model dependence. 
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Figure 1: Mg versus V: The pink region comes from the constraint of the life time of the g-quark, which 
should be less than 0.1 sec. The green region comes from the constraint of the proton stability. The black 
region is allowed by the both constraints. The heavier of Mg, the wider the allowed region is. 



4 The Maki-Nakagawa-Sakata matrix 

In this section we confirm the assignment of Table 2 realizes the neutrino masses and MNS matrix. In the 
superpotcntial 

Wl = Yf [H^iLiN^ + L2NI) + H^{LiNl - LsTVs')] 

+ Y^El{H^Li + iff L2) + Y^E^H-^Ls + Kf £;^(i7f L2 - H^Li) 
+ ^Y,^'^3iN!N^ + N^N^) + ^Yi'<P3N^N^ 

+ iy/^[2<i>i7Vf7V2^ + $2(7Vi=7V^=-iV2^7V=)], (27) 
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we define the VEVs of scalar fields as follows 

($1) = |a|-l^oCAr, ($2) = |a|K)Sw, ($3> = Vb = ^^=, (28) 

VI + lap 

and define the mass parameters as follows [14] 

~-Y,^Wo, M2 = Y^^''\a\Vo 

(29) 

Without loss of generality, by the field redefinition, we can define Y^2 31 ^1*3 j ^^^-l ^^'^ non-negative 

and Y2 '^ , Y^ are complex. For simplicity, we put 

yM^yM^-^^ yM^giV^ (30) 

and 

a = e"''|a|, Mi = M3 = V^,, M2 - 0^0- (31) 
With these parameters, the mass matrices are given by 

Ml = 4^\ m\ j , Md^^ 

1 + asjv clcn 
Mr = Vq ( acN 1 — asjv 
1 

Due to the seesaw mechanism, the neutrino mass matrix is given by 
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The charged lepton mass matrix is diagonalized as follows 

V^MlMlVi = diag{ml,m'.,,m',)^{{m'2)\{rn\)\{m\f), (34) 



Vi ^ —\ 1 1 . (35) 



To realize experimental results [15], the neutrino mass matrix should be diagonalized as follows 
Vl My Vy = diag{my^ , m^^ , m^^) , 

(-CyCiz -SyCis Si3 

-75"- + 7I^-"l3 TfC, + -^55,513 ^Ci3 

sin^ 20^ = 0.8704, sin^ 26I13 0.11 

-ml^ - -ml, = 7.6 x 10"^ [eV% ml^ - ml^ ^ 2.4 x 10"^ [eV^]. (36) 

Due to the overabundance of parameters, unfortunately, it is impossible to fix the parameters by these 
constraints. Therefore we assume a and p| are real for simplicity, then we get 

a = 1.40321, SAT = 0.110194, 

= 0.0212781 [eV], = -0.0539146 [eV], P4 = 0.112142 [eV], (37) 
TO^i = -0.0207554 [eV], m^^ = 0.0225119 [eV], m^, = -0.0539146 [eV], 
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and 



AY/^Vu)'^ _ /I c/i 1nl2^ fvN 



y=v/rTH^^^i-J^ = (1.54xl0i^)(y/^^j [GeV]. (38) 

From the requirement of pcrturbativity of Yukawa coupling such as Y/^ < 1, y is bounded from above. It is 
difficult for V to be much larger than O(lO^^GeV). Therefore, there exists upper bound for V even in the 
case that proton decay is perfectly forbidden. It must be noted that there is non-trivial coincidence of the 
constraint of Eq.(26) with RHN mass scale. 

Finally we give a comment about how to realize the VEV directions given in Eq.(28). The Higgs potential 
derived from superpotential 

+ X^H^iSiH^ + S2H^) + X^,{SiHY + S2H^)H^, (39) 

has accidental 0{2) symmetry. To avoid massless Nambu-Goldstone boson, we add soft S4 x Z2 breaking 
terms in the form of the inner products with (1,1) as follows 

£ D ml^{Hl,')\H^ + H^) + mU{Hi)\H? + H^) + m|s(^3)^(5i + ^2) + /i.e., (40) 

then the VEV direction (^1,^2) oc (1,1) {A = S,H^ ,H°) becomes the minimum of potential and the 
VEV direction of , , Sa in Eq. (28) is realized. The potential of $q, $^ derived from the leading order 
superpotential 
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2Mp 



{2^/3$l$2 [mf - m'] + i'^i - *2) [m? + - 2{<^i?] } 

$3 {V3<^i [mf - mf] + *2 [mr + ~ 2(a>D'] } , (4i) 



does not have accidental symmetry. To avoid domain wall problem, we must add soft S4 breaking terms. 
Then the VEV direction of $, is controlled by the parameters 2,3.4 and soft SUSY and flavor breaking 
parameters. The mechanism for inducing soft flavor symmetry breaking terms is unknown and beyond the 
scope of this paper. We leave it for future work. 

5 The origin of the scale of V 

Finally we explain how the required value for V is realized. The superpotential Eq.(41) is simplified as 
follows 

= ]^(**')'- (42) 

The origin of ^-potential becomes unstable point due to the negative soft SUSY breaking squared mass and 
the potential is lifted by F-term derived by Vlinimizing the potential, the VEV of $ is estimated as 

follows 

(43) 

For the typical range of A and SUSY breaking scale msusY such as 0.01 < A < 1, 0.1 TcV < msusY < 
10 TeV. Hence the region of V is as follows 

10^° GeV < < 10^2 Qgv. (44) 

Although the each of region given in Eq.(26) and Eq.(44) is very narrow, remarkably, there exists overlap. 
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6 Conclusion 



We have considered the suppression mechanism of proton decay based on 6*4 x Z2 flavor symmetric model. 
Under the field assignment that MNS matrix is realized, we have classified the several suppression mecha- 
nisms. There are two new solutions other than the well known leptoquark solution. For the cancellation 
solution, the four-Fermi interaction which induces proton decay vanishes in a special VEV direction. For the 
suppression solution, the stability of proton is satisfied for appropriate size of VEV. Although the allowed 
region for the VEV is very narrow, there is coincidence between the allowed regions required by the differ- 
ent phenomenological considerations such as naive potential analysis, RHN mass scale and the life times of 
g-quark and proton. 

A Mixing matrix of scalar g-quarks 

Here we define the mixing matrix of scalar g-quarks. The mass terms of scalar g-quarks are given as follows 

-C D mU\Gi\' + \G2\' + \Gr)+^^U\G'i\' + \G^2\' + \Gl\') 

+ kAk[S3{GiG1 + G2G'2 + G3G^) + h.c] 

+ |fc(GiGJ + G2G2 + GsG's) + XiiH^H^ + H^H^) + X^H^ H^\^ 

+ \k\^\S3\\\G,\^ + \G2? + |G3|2 + |GJ|2 + |G^|2 + \Gl\^) + D-terms 

= I](G:,G^) ^ ^^^^y ^ , (45) 

where we assumed Ak is real for simplicity. If 

Ml_ = kAkv'^ + k{XiVuVd + Xsv'y^) = 0, (46) 

is satisfied, then there is no G — G'^ mixing. In the case that 0, the mixing matrix of scalar g-quarks 

is defined as follows 

Vo^{ Z 7; ) , ( ^ ) V, ^ diag(Mi,M5). ( )^Va{ ) , (47) 

From this definition, the propagators of scalar g-quarks are given by 

where Mq = a/M+MI is mean scalar g-quark mass. 
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